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ABSTRACT 

Recent studies by iLada et all (|2010f ) and iHeiderman et al.l ()2010f ) have suggested that efficient star 
formation occurs above an approximate threshold in gas surface density S of Sc 120 Mq pc~^ 
{Ak ^ 0.8). We find no precise threshold for star formation; the impression of such results from a 
continuous and steep power-law increase of the ratio of protostellar mass to molecular gas mass with E, 
approaching unity at protostellar core densities, corresponding to S ~ 1000 M© pc~'^. We argue that 
this increase in star formation efficiency results from the increasing importance of self-gravity with 
increasing density, along with the consequent decrease in evolutionary timescales. The observations 
are consistent with models in which regions of more diffuse molecular gas with column densities 
corresponding to Ay ~ 1 — 2 are initially formed by converging galactic hydrodynamic flows which 
subsequently collapse gravitationally, producing a power-law relation between surface density and the 
area A spanned at that density of ^4 ~ We show that the finding of a strong correlation between 

the amount of gas above Ec and the young stellar population by Lada ct al. (2010) with ratios q of 
dense gas mass to stellar mass of 8 > q > 3 requires continuing formation of dense gas, also consistent 
with the idea of gravitational collapse. Finally, we note that the suggested linear relationship between 
the star formation rate and molecular gas mass at very high surface densities in galactic studies arises 
as a result of probing small size and high density scales where rapid gravitational collapse is occurring 
and the efficiency of star formation approaches unity. 

Subject headings: ISM: clouds, ISM: structure, stars: formation 



L INTRODUCTION 

Local star formation has long been known to be 
strongly en hanced in the densest region s of molecular 



clouds fe.g.. iLadal 



[19951 lOnishi et all 



1992 



199a 



Lada et al.lll993; iMizuno et al.l 
Johnstone et al.ll2004f ). The in- 



creasingly large and sensitive surveys of molecular clouds 
using a variety of techniques have now made it possible 
to derive more quantitative relationships between star 
formation and dense gas. I n particular, using measure- 
ments of infrared extinc tion, ILada et al.l ()2010l LIO) and 
IHeiderman et al.l ()2010l . HIO) proposed that galactic star 
formation is strongly enhanced above a critical gas sur- 
face density corresponding to S,. ^ 120 M^ pc~^, ap- 
proximately the same threshold lOnishi et al.l ()1998[ ) de- 
rived from C^^O observatio ns of the Taurus molec ular 
cloud. In addition, HIO and ILada et~aI1 ()2010l [20121 ) in- 
ferred that star formation rates (SFRs) depend approx- 
imately linearly on the amount of molecular gas above 
the threshold, qualita tively consistent with findings from 
extragalactic studies (|Gao fc Solomonir2004l : iBigiel et"all 
1201 Ih . although with differing normalizations. 

In this paper we suggest that there is no precise thresh- 
old for star formation, but densities must be high enough 
so that the self-gravity of the gas is important. We 
show that the observations are consistent with a contin- 
uous progression of the increasing importance of gravity 
and decreasing evolutionary timescales with increasing 



density. We argue that the observations of LIO gener- 
ally require continuing formation of dense gas, as pre- 
dicted by models of molecular clouds undergoing large- 
scale gravitational collapse (e.g., Har tmann et al., ,2001 
Vazqucz-Scmadcni ct al 
Heitsch 



2007 
Hartmann .2008. 
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Zamora-Aviles etal\ l2012f) . 
We point out that the relation between the star forma- 
tion rate (SFR) per unit area and the surface density 
found by HIO only becomes linear at very high surface 
densities where the efficiency of star formation is very 
large; we suggest that this results in a linear relation be- 
tween the star formation rate and S, particularly when 
density-sensitive tracers are used. Finally, we show that 
the data of HIO suggest that at high densities clouds ex- 
hibit a similar power-law dependence of the area A at a 
given surface density E, A{T,) oc which might pro- 

vide important constraints for theories of the structure 
of molecular clouds. 

2. OBSERVATIONAL EVIDENCE FOR THRESHOLDS 

It is useful to first review the observational constraints 
in some detail. We begin by considering the findings of 
LIO who examined the ratio of the gas mass above a given 
level of extinction (or surface density) Mg to the number 
of young stars N in several nearby star-forming regions. 
They searched for the value of E which would most nearly 
result in a linear relation between Mg and N for the 
various clouds; the resulting critical value was found to 
be Ec 116 ± 25 M0 pc'^. While LIO identified this as 
a threshold for star formation, they did note that there 
was a factor of nearly two in E for which the dispersion 
in Mg/N among the various clouds was minimal. 

To examine further the evidence for a threshold surface 
density, we consider the findings of HIO. These authors 
established roughly equally-spaced contour levels of (in- 
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frared) extinction in a sample of low-mass star-forming 
regions, and then measured both the area and the num- 
ber of protostellar (Class I and flat-spectrum) sources 
contained within each extinction or surface density con- 
tour. By focusing upon objects thought to be protostars, 
and thus a) very young and b) still accreting from their 
natal envelopes, HIO tried to minimize the conceptual 
difficulties which arise because the gas one observes at 
the present epoch is not the gas that produced the stars. 

The left panel of Figure [1] shows the main results of HIO 
for the low-mass star-forming regions. We plot this as the 
number of young stars N per unit contour area A vs. S; 
this is equivalent to the star formation rate of HIO, which 
is simply N/A multiplied by an average protostellar mass 
and divided by a (constant) timescale of 2 Myr. (We have 
added together the Class I and flat-spectrum sources to 
improve the small number statistics; see discussion in the 
Appendix.) The ratio N/A is steeply dependent upon the 
molecular gas surface density, increasing by about three 
orders of magnitude over a range of about one order of 
magnitude in E, qualitatively consistent with the idea of 
a threshold, although there is no sign of a transition in 
the N/A vs. S relation at their suggested threshold of 
Ec = 129 ± 14M0 pc~2 (dotted vertical line in Figure [l]). 

In the Appendix we discuss the findings of HIO in 
more detail. Here we note that the correlation between 
N/A and E is mostly driven by the correlation between 
the area A of the contours as a function of E, where 
A oc E~^'^. Thus, the order of magnitude range in E 
probed in Figure [1] corresponds to a range of three orders 
of magnitude in A; moving to higher values of surface 
density implies restricting attention to ever-decreasing 
areas and volumes. In addition, if N is independent of 
E, N/A increases by three orders of magnitudes, consis- 
tent with the trend, shown in the left panel of Figure 

In fact, the relevant physical quantity for star forma- 
tion is more likely to be the volume density p than E. If 
we make the guess, as discussed further in the Appendix, 
that the typical length scale along the line of sight at E 
is I ~ A^/^, the range of one order of magnitude in E 
spanned in Figure [T] corresponds to a range of roughly 
a factor of 300 in p and a decrease by a factor of 17 in 
collapse timescale. It is therefore not surprising that the 
number of young stars per unit area, or the star forma- 
tion rate per unit area as in HIO, is a rapidly increasing 
function of E (see discussion in [J3]). 

The impression of a threshold for star formation is en- 
hanced in HIO by their inclusion of an additional data set 
spanning higher surface densities than present in their 
study of low-mass regions. Including data from mas- 
sive clumps with surface densities measured from HCN 
and star formation r a.tes SFR estima ted from infrared 
luminosities Ljr from lWu et all ()2010[ ). HID find a more 
nearly linear relationship between the SFR and E, span- 
ning the range from 300 < E/Mgpc-^ < 3000. The 
intersection of this linear branch with the low-mass star 
formation branch, discussed above and shown in Figure 
[U lies at about Ec = I2OM0 pc~^ which HIO interpret as 
a signature of a threshold. However, as HIO acknowledge, 
there are concerns about combining results using two dif- 
ferent methodologies, from two very different types of 
regions (low-mass star-forming regions vs. massive cores 



with Ljfi > 1O'*'^L0 and thus forming massive stars) (see 
also fj5]). First of all, the physics of high-mass star forma- 
tion might differ from low-mass star formation. Secondly, 
the results from the low-mass regions do not identify a 
particular value of Ec but instead extend smoothly into 
regions with E > Ec. As a resul, the linear branch of 
high-mass star forming clumps does not start where the 
low-mass star formation branch ends. It instead inter- 
sects the low-mass branch somewhere as one would ex- 
pect if both branches are two independent correlations 
with different slopes. In this case, the surface density Ec 
at the point of intersection has no physical meaning and 
cannot be interpreted as a threshold above which star 
formation depends linearly on gas surface density with a 
steeper dependence below this threshold. 

Assuming for the moment that high-mass star forma- 
tion proceeds in a similar way as low-mass star forma- 
tion, how can we then reconcile the linear behavior of 
the high-mass star formation branch with the strongly 
non-linear relationship, observed for the low-mass star 
forming regions? The linear, high-mass star formation 
branch traces gas with molecular core densities of > 10^ 
cm^^. To examine the implications for low- mass star for- 
mation we examine the star-gas correlation at these high 
densities. In the right panel of Figure [T] we show the ra- 
tio of protostellar masses to the mass of gas within the 
contour at E 0. At high surface densities the "efficiency" 
of star formation, as measured by the ratio of dense gas 
mass to protostellar mass, exceeds 0.1 and reaches val- 
ues of order unity at molecular core surface densities of 
E ^ 1000 M© pc^^. At this point the non-linear increase 
of stellar mass versus dense gas mass with surface den- 
sity is expected to break down as it is unlikely that more 
than 50% of a core turns into stars. If there exists a max- 
imum or typical efficiency for dense cores to form stars 
one would therefore naturally expect to observe a transi- 
tion towards a linear N/A vs. E relation above a critical 
value of Mt:/Mg. As suggested in the Appendix, if in the 
massive cores the HCN tends to trace these particular, 
dense environments and thus also a particular dynamical 
collapse time, a linear relation between the gas mass and 
star formation rate might ensue. 

One should also note that the area spanned at the high- 
est observed surface densities is of order A ~ 10~^pc^ 
(see Appendix). This implies that the length scales 
probed at the highest values of E are ^ A^^^ ~ 0.3 pc. 
The HIO data is averaged over 270 arcsec beams, which 
at the typical distances of the studied clouds of 150 — 
300 pc corresponds to resolutions of 0.2 — 0.4 pc. Thus 
dense, low-mass protostellar cores with typical sizes of 
0.1 pc (jMvers fc Bensonlll983 ^. wiU drop out of the high- 
E region; cores that are surrounded by much lower den- 
sity material will then appear at lower- E. This again 
emphasizes the conceptual difficulties involved in relating 
the (extended) gas seen at the current epoch to the gas 
that formed (or is most directly forming) the observed 
stars. 

3. "THRESHOLDS" VS. CONTINUOUS EVOLUTION 
3.1. Overview 

* We caution that this plot shows intrinsically correlated pa- 
rameters, as the gas mass is derived by multiplying the area of the 
contour by S, so the slope of the relation is not significant. 
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As shown in the previous section, while HIO find a 
rapid increase in the SFR or N/A over a region in sur- 
face densities that includes the suggested threshold sur- 
face density, there is no obvious break in the relation, pre- 
cisely at Sc- Still the rapid increase of M^/Mg for surface 
densities that lie within a factor of 3-4 around Ec could 
be interpreted as a threshold for efficient star formation. 
Consistent with these results, here we argue that appar- 
ent thresholds like this occur in continuously-evolving 
clouds because: most molecular clouds form at modest 
densities and pressures; the gas must become much more 
dense on its way to form a star; the evolutionary (gravita- 
tional collapse) timescales rapidly become shorter as the 
gas becomes denser; and therefore the contrast of low- 
density, slow evolution vs. high-density, fast evolution 
leads to the impression of a threshold. 

Consider an estimate of the central pressure of a gas 
cloud P = Pex + ttGT? /2, where Pex is the external 
gas pressure and S is the su rface or column density 
(jElmegreen fc Elmegreenlll978f) . This relation holds ex- 
actly for an infinite hydrostatic sheet and is roughly 
correct for many other geometries. Now, at the sug- 
gested critical surface density of Sc ~ 120Mopc~^, 
P ^ 100 X Pel if we assume a typical local interstellar 
medium (ISM) pressure Pex/k ~ lO** cm~^ K. 

To achieve such ram pressures from flows at lOkms"^, 
typical of the turbulent I SM velo city dispersions in regu- 
lar star- forming galaxies (iDib et a l. 2006) , would require 
very large external flow densities of > 35cm~^. In addi- 
tion, to produce a region of this extinction, correspond- 
ing to hydrogen column density o{ Nh ^ 1.4 x lO^^cm"^, 
would require the converging or expanding flows to last 
for a timescale 

S , , /I0kms-i\ /2cni-3\ 
r = = 1.1 X lOV , 

2Woo?^oo V Voa J \ Tloa J 

which is a factor of ten to twenty times larger 
than expected from large-scale numerical simulations 
of t he dynamics of th e turbulent ISM in disk galax- 
ies (jPobbs et al.l I2012D that indicate cloud formation 
timescales of order 5 x 10^ — 10^ yrs; again, one would 
require initial densities of 20 — 40 cm~'^ to build the cloud 
in reasonable timescales. 

Most of the mass of molecular (CO) clouds in the so- 
lar neighborhood lies at column densities co rresponding 
to Ay ~ 1 - 2 (e.g., [Goldsmith et al.ll2008f ): i.e., there 
is an extended, lower-density molecular envelope sur- 
rounding the dense, strongly self-gravitating gas, which 
would be at best slowly contracting due to its self grav- 
ity (see below). At Ay ^ 2, or a surface density 
~ 4OM0 pc~^, the pressure due to self gravity would be 
P/k >6x lO^'cm"'^ K, only about an order of magnitude 
larger than the typical ISM pressure. This is indeed a 
natural result of forming molecular clouds (at least in 
the solar neighborhood) by interstellar medium (ISM) 
fiows with densities of > a fe w cm~^ at v elocities of or- 
der 10 km s~ ^ , as simulated by (jBall ^tCTO S-Paredes et alj 
1999; Audit & HenncbcUc 2005; Hcnncbellc & Audit 
2007i: iVazguez-Semadeni et al.ii2007i : iHeitsch et al..,2008t : 
Heitsch fc Hartmannll2008l) . Moreover, such surface den- 
sities can be formed in the requisite timescales, especially 
if the inflowing material has slightly higher densities than 



that of the average ISM ()Dobbs et al.ll2012[ ). 

Thus, while the majority of gas in molecular clouds 
is at (relatively) low densities, and can be produced 
by external ISM ram pressures, gas at or above Ec 
is very likely to be generated and bound mainly 
by gravity. Independent support for this argument 
comes from analyses of the probability density func- 
tions (Npdf; the equivalent of the A(E) relations dis- 
cussed in §2 and in the Ap pendix) of the column 
densi ties of molecular clouds (" Kainulainen et all I2009L 
120111 : Trocbrich & Rowlcs 2010). The Npdfs appear 
lognormal at low densities, which is generally inter- 
preted as the res ult of turbulence-dominated flows 
(|Vazquez-Semadeni|[l99llOstriker et al.ll200lD . but they 
exhibit power-law tails at high column densities in star- 
forming clouds, as do the result s of H IO (§2, Ap- 
pendix). 'Ballesteros-P aredes et al.l (|2011[ ) showed that 
such power- law tails naturally arise when gravitationally- 
driven motion becomes more important than pure hydro- 
dynamic turbulence, consistent with our picture. 

Note that although HIO find a rapid increase in star 
formation efficiency of about three orders of magnitude 
near Ec, it actually takes place over a range of one or- 
der of magnitude in E which corresponds to two orders 
of magnitude in "self-gravity pressure" . Moreover, we 
argued previously that this range in E corresponds to a 
factor of ~ 300 in volume density; this corresponds to 
a decrease in the free-fall time oc p^^l"^ of more than 
one order of magnitude at high densities. Thus the com- 
bination of strong self-gravity and rapid evolution will 
enhance the impression of a threshold. 

3.2. Self- gravitating sheet simulation 

To illustrate our argument in a more graphical way, 
we use some simple numerical simulations. We con- 
sider a massive sheet of gas which could have arisen 
as a result of colliding fl ows a s in the simula t ions o f 
lyazQuez-Semadeni et al.' C2007I); iHeitsch et aD (|2008| ): 
Heitsch & Hartmann (2008). We ignore turbulent mo- 
tions and assume that the cloud is initially completely 
quiescent but generates supersonic velocities via global 
and local gravitational co llapse. The setup is basicall y 
the same as that used in iBurkert fc HartmannI (|2004l ). 
where we considered a uniform circular sheet. 

The sheet is in hydrostatic equilibrium and in pres- 
sure equilibrium with a constant surrounding pressure of 
P/k = 5 X 10^ K cm~'^. The initial conditions are the ra- 
dius R of the sheet, its mass M, and the sound speed Cg of 
the gas. We assume an isothermal equation of state with 
the pressure P = c^p. Here, we took an initial cloud ra- 
dius of i? = 10 pc, a sound speed of Cg — 0.2kms~^ and 
calculated results for initial sheet masses ranging from 
Mc = IO^Mq to 2 X IO^Mq. The calcul ations were per- 
forme d wit h the SPH code o utline d by iBate fc BurkertI 
(5997i) and IBurkert &: AlvesI ((20091) . 

Figure [2] shows the surface density of our standard 
model with Mc — at two times. With the adopted 

initial conditions, a dense ring develops near the edge 
of the cloud, due to highly non-linear a ccelerations of 
gravity fe.g.. IBurkert fc Hartmam] 120041) . In effect, we 
have made a circular and in this case artificially smooth, 
dense filament, situated within a lower-density cloud. 
The ring/filament continues to grow as the cloud (which 
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has many Jeans masses) globally contracts, accumulating 
mass and becoming denser as time proceeds. 

Figure |3] (left) shows the evolution of the mean ring 
volume densities for the standard model (solid line) as 
well as for additional models with differing initial masses 
and thus surface densities. The density evolution can 
be characterized by two phases, "slow" and "fast" , a dis- 
tinction which is clearest for the two lowest-mass models. 
The upper limit to the "slow" ti mescale is the collapse 
time tcoH for the sheet as a whole (jBurkert fc HartmannI 
12004,), 

1/2 

« I ^ I ^ (G(p))-i/' . (2) 



In other words, tcoii is the global free-fall time of the 
sheet characterized by (p), the mean density over the 
spherical volume enclosed by R. Equation ([2]) yields a 
global collapse time of 15 Myr for Mc = IO^Mq. In 
this case the sheet contracts considerably before the ring 
(filament) density runs away. Because the sheet volume 
density scales as , the timescale for collapse of the ring 
becomes much faster with increasing surface density. 

Figure [3] demonstrates our basic explanation of ap- 
parent thresholds for star formation. The timescales 
over which the cloud remains at low column densities 
are one to a few Myr, consistent with the global sheet 
collapse timescale. The surface density in the ring at 
first grows slowly as it is initially pressure supported. At 
some point however enough material has been accumu- 
lated and the ring becomes gravitationally unstable. It 
collapses onto itself, leading to runaway growth in den- 
sity that is much faster than the growth in mean density 
in the sheet and resulting in two gas phases, a diffuse 
gaseous component and an embedded dense collapsing 
filament that now can fragment and form stars. The col- 
lapse occurs essentially at the free-fall time of the local 
density, very similar to the evolution under pure free-fall 
that is shown in the right panel of Figure [31 The two 
timescales of the collapsing ring and of the diffuse sheet 
are comparable to observational estimates of cloud life- 
times (Hartmann ct al. 2001) and protostcllar/starless 
core lifetimes (Eva ns et al.ll2009i ). respectively. 

Of course these simulations are also highly unrealistic, 
in that clouds will have multiple sites of star formation; 
but one might view Figures [2] and [3] as representing the 
evolution of differing parts of the cloud or even differ- 
ing local patches, with the same basic result. Indeed, 
even the limiting case of the uniform pressureless sphere 
collapse exhibits the aspect of slow initial growth and 
localised fast runaway collapse. 

Note that in our simulations, the entire cloud is grav- 
itationally bound and collapsing. However, even if we 
do not support the low-density medium via e.g. tur- 
bulence or magnetic fields, a broad range of evolution- 
ary timescales d evelops due to ring formation (see also 
IPon et al.l[20Tll) . 

4. EFFICIENCY OF STAR FORMATION ABOVE Sc 

LIO argue that the total number of young stars, or the 
stellar mass , is much better correlated with the mass 
of dense gass Mg above the threshold Ec than with the 
total gas mass in the star-forming region. To explore 
this further, in Figure |4] we plot the data from Table 2 



of LIO in the form of the ratio of the dense gas mass to 
stellar mass Mg/M* = q, vs. the ratio of the total gas 
mass Mtot to Mg, again using an average mass per star of 
O.7M0. Discounting the Lupus 4 region for poor stellar 
counting statistics, and noting that the stellar population 
of Orion B may be underestimated due to high extinc- 
tion and nebulosity (S.T. Megeath, personal communica- 
tion), the data suggest that the least- active regions with 
large values of q have low fractions of dense gas, while in 
the most active regions the fraction of dense gas is high, 
which may be an indication of cloud evolution. For the 
most active and probably evolved star-forming regions, 
(apart from Orion B) there is a relatively modest range 
in the gas to stellar mass of 8 > g > 3. 

An important implication can be derived from these 
results. As all the young stars have not formed at the 
same time (indeed, some of them are still-accreting pro- 
tostars) , must be increasing with time; therefore the 
mass of dense gas Mg must also increase in time to keep 
q above unity and within a modest range. Thus, the LIO 
results imply that star-forming clouds are n ot in a quasi- 
static state but are continuing to evolve (jPobbs et al.l 
2011a,b, 20W)- 

To illustrate further implications of the data, we con- 
sider a simple analytic model. Let us assume that a 
dense molecular gas component is generated at a rate 
Min- Suppose also that a fraction e of the dense gas 
turns into stars on its local free-fall timescale r/ / 



dMg dM, 



dM^ 
dt 



-M„ 



(3) 



(4) 



where we take e and Tff to be constants. Typical values 
of e are e « 0.1 (LIO). Let us first demand that at any 
given time t the ratio q = Mg/M^, is constant. We then 
find 

Mg{t) ^ Mg,Q expit/to) , (5) 

where Mg^ — Mg{t — 0) is the initial dense gas mass 
when star formation starts and io is the exponential 
growth timescale of the dense gas mass. 



^0 = ^. 

Inserting equation ([S]) into equation (jS]) we find 



M,„ =e(g + l)- 



M, 



9,0 



Tff 



exp 



(6) 



(7) 



Thus, if q were absolutely constant, there would need to 
be an exponential infall and increase of dense gas mass 
with time. 

This exponential growth model raises some concerns. 
First of all, obviously Mg cannot increase exponentially 
forever. In addition, as the fraction of dense gas is ob- 
served to be always small compared to the diffuse gas 
mass and as the efficiency of star formation is gener- 
ally only a few percent of the total clou d mass (e.g., 
lEvans et al.l l2009t iKrumholz et al. I I2012D , exponential 
growth would imply a very delicate timing for disper- 
sal of the whole cloud. Stellar feedback would have 
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to destroy the cloud before a large fraction of its dif- 
fuse gas has been converted to dense gas. Finally, it 
is not clear that simulations provide support for con- 
tinued exponential growth of the dense gas mass (e.g., 
I Vazquez- Semadeni et al ."2009'. see §3.2|) . To address this 
problem we develop a second analytic model which as- 
sumes linear rather than exponential growth in the dense 
gas mass with time. Adopting time units r = t/rff and 
assuming a constant Mm, the solutions of equations ^ 
and (m now are 

Mg{T) = Mg,o e-^" + ihn T//(l - e-^")/e, (8) 

M,(t) =Mg,o(T-0)(l-e-^-) 

+ M„.r//[r - (1-e— )/e]. (9) 

Figure [5] shows results, assuming Afg.o = for efficien- 
cies of 0.1 and 0.3, respectively. For values of e w 0.1 that 
are consistent with the observations of TIO, q remains in 
the observed range over 4t//. 

LIO suggested that the critical surface density might 
correspond to a critical volume density ric = 10^ cm~^, 
implying Tff ^ 0.35 Myr. Thus, with linear growth of 
the dense gas mass. Figure [5] indicates that even for star 
formation efficiencies of e = 0.1, the observed range of q 
can only be maintained for less than 1-1.5 Myr, consider- 
ably shorter than typical estimates of star-forming cloud 
lifetimes which are of order 2 — 4 x 10^ yrs. However, 
there are reasons to question this estimate of the criti- 
cal volume density. In the Appendix we use the data of 
HIO to examine the issue further, by estimating the path 
length along the line of sight to be I ^ A^l"^ and thus 
< p >^ E/A^/^. While the results differ from cloud to 
cloud, we find a narrow range of estimated mean densities 
at Ec between ric ~ 500cm~^ and ~ 2.6 x lO^cm"^. Tak- 
ing a mean value for all clouds of ric ^ 10^ cm~^ implies 
a free-fall time of 1 Myr. This makes it much easier to 
explain the range of q found by LIO as q now would lie in 
the observed regime for 4 x 10^ yrs which are resonable 
star formation timescales of molecular clouds. As a con- 
sistency check, the areas at Ec of one to several pc^ are 
much larger than those typical of low mass pre-stellar 
cores of densities of 10^ cm~'^ (e.g., iMvers fc BensorJ 
[1983). 

We can also use the numerical simulations to check the 
arguments constraining q. The left-hand panel in Figure 
[6]shows the growth of mass with densities above lO'^cm"'^ 
for the standard model. One observes that the growth 
of dense gas is exponential for the first ^0.2 Myr, but 
then becomes roughly linear. If we were to keep q strictly 
constant, the cloud would have to be dispersed at the end 
of the short epoch of exponential growth which is not in 
agreement with observed star formation timescales. To 
illustrate the implications further, in the right panel of 
Figure m we show q = Mg/M^, as a function of time for 
the standard model with the mass in the circular 
filament above n{H2) — lO'^cm"'^, but at a free-fall time 
of 1 Myr earlier, times an efficiency factor e. In other 
words, here we assume again that tMg turns into stars 
after one free fall timescale. The dotted curve assumes 
e = 1, the solid red curve e = 0.3, and the dashed curve 
e = 0.1. 

Figure |6] reinforces the analytic results. For efficien- 



cies of order 10% q will lie in the range found by LIO 
for 4 free fall timescales. If, in addition, the free fall 
timescale is of order 1 Myr, as expected for critical den- 
sities of n{H2) = lO'^cm"^, q would remain in the ob- 
served regime for the typical lifetimes of star-forming re- 
gions. The pure analytic model actually does a good job 
of reproducing these combined numerical- "post process- 
ing" results in Figure [6] (for example, it predicts that for 
e = 0.1, q drops below ~ 3 at Arff ~ 4 Myr). 

Dense gas formation is probably not monolithic, but 
instead is the result of differing regions collapsing at dif- 
fering times. The blue solid curve, labeled /, shows the 
situation if we assume that after every t// = IMyr a 
new, independent dense gas region forms that does ex- 
actly the same thing as the e = 0.3 case, just starting 
later. So e.g. after 2 Myr we have in total 3 star forming 
regions: one is 2 Myr old, one is 1 Myr old and just be- 
ginning to form stars, and one is just starting to generate 
dense gas but is not forming stars yet. This spreading of 
the onset of star formation is helpful but does not change 
the main result. 

Finally, the behavior of q implies something about 
the efficiency of turning gas into stars. Sometimes star 
formation sites are characterized by the "star forma- 
tion efficiency per free- fall time" rjff = SFR x Tff/Mg 
(|Krumholz et al.1l20l2[ ). Considering only the gas mass 
above Ec, LIO find a low value r/ff « 0.02; however, 
this depends upon their assumption of ric = lO^cm^**. If 
instead ric = lO^cm""*, then r]ff ^ 0.06. 

The concept of ?7// is complicated, not only because 
real clouds exhibit a wide range of densities and thus re- 
gions with very different free-fall timescales, but also be- 
cause it depends upon the assumption of a quasi-steady 
state. For example, if the gas mass grows exponentially 
as in the model which provides strictly constant q, 

vff = (qt/Tffy^ . (10) 

Note that rjff, despite its name, does not depend ex- 
plicitely on the star formation efficiency e, that is the 
fraction of dense gas that turns into stars in r^j . A low 
value oi rjff is therefore not in confiict with theoreti- 
cal models that would predict higher values of e. Given 
q, rjff is instead a measure of the lifetime t of the star 
forming region in units of its free fall time that is being 
continuously fed by infall from its diffuse molecular enve- 
lope. Even if the growth is not strictly exponential, there 
is still a tendency to underestimate r]ff, as well as signif- 
icant uncertainty as to what t he conditions were in the 
past when the first stars formed lVazauez-Semadeni et al.l 
(e.g.,[200i. 

5. DISCUSSION AND CONCLUSIONS 

We have shown that the observations of HIO and LIO 
do not require precise density thresholds for star forma- 
tion; all that is required is rapid gravitational collapse at 
high densities coupled with the presence of a much lower- 
density molecular cloud formed via plausible interstellar 
medium flows. This "external" low-density molecular 
material feeds the dense regions. The low-density gas 
can also be gravitationally collapsing and still provide a 
more slowly-evolving structure with more mass than in 
the dense regions, especially if the global cloud geometry 
is far from spherical. If the flows into the dense regions 
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are driven by gravitational acceleration, the increase of 
dense gas mass with time to keep the ratio of dense gas 
to stars relatively constant naturally occurs. The result- 
ing picture is consistent with simulations of dynamically- 
evolving, gravitationally-coUapsing star-forming molecu- 
lar clouds. Density diagnostics are needed to help trans- 
late surface densities into volume densities, which then 
will constrain the free-fall times and thus the efHciences 
of star formation at modest surface densities. 

The data of HIO suggest that, at high E, the efficiency 
of converting the mass of gas into stars appears to ap- 
proach unity; we suggest that the linear relationship be- 
tween SFR and S then arises because this dense gas is 
already in gravitational collapse. The survival of this 
linear relationship in extragalactic, beam-diluted obser- 
vations suggests an approximate common scaling of the 
mass of dense gas with gas mass above Ec; an open 
question is whether this is a result, or can be derived 
from, the apparent "power-law" distribution of areas as 
a function of density. The origin of the A{Y^) relation 
is unclear and has to be related to the universal den- 
sity structure of clouds; the transition between the log- 
normal form of the Npdf at low column densities to the 
power law behavior seen in studies of star formation by 
HIO a nd others is likely the product o f gravitational col- 
lapse (|Ballesteros-Paredes et al.ll2011[ ). 

Along the way we have noted some of the conceptual 
difficulties in deriving density thresholds for and efficien- 
cies of star formation. The dense gas we see at the 
present epoch is not the gas that formed, or is forming 
the young stars currently present. Can we then assume 
a steady state to interpret the observations and linking 
past star formation to present cloud structure? In the 
case of galactic molecular clouds, we have argued that 
one in general cannot assume such a steady state. Be- 
yond this, will the gas that we see at any density at the 
current epoch turn into stars later? Yet another con- 
ceptual problem is what volume(s) should one choose to 
examine in order to determine a threshold for or effi- 
ciency of making stars. In the picture we have presented 
here, with continuous evolution of lower-density gas into 
high-density, star-forming gas, it is not obvious how one 
would choose any particular scale. As stars can only 
be formed from gas that has much higher densities than 
that of the interstellar medium in general, it is not clear 
what one learns from connecting young stars with the 
very densest gas. Extragalactic observations over large 
scales may sidestep these questions to some extent by 
averaging over sufficient volumes that steady states may 
be achieved; but then the question is whether physical 
insight is suppressed by such global averaging. A fuller 
understanding of star formation must involve the cycling 
of gas between low- and high-density phases. 

This work was supported in part by NSF grant AST- 
0807305. AB acknowledges useful conversations with 
Mark Krumholz and Clare Dobbs, LH with Javier 
Ballesteros-Paredes. 

6. APPENDIX 

Here we look a little more closely at the data which 
go into producing the HID results illustrated in Figure [U 
using the data from their Tables 1 and 2. 



While the N/A vs. E relation appears well-defined, it 
is in fact more of a reflection of the behavior of A vs. E 
than of the protostars. This can be seen in the left panel 
of Figure [7] which shows a weak correlation between the 
number of protostars in a contour at a given E. The 
correlation shown in Figure [T] is driven by the area in a 
specific contour A vs. E behavior, which is shown in the 
right panel of Figure [T] Of course protostars are formed 
necessarily by high-density gas; the lack of correlation in 
the left panel is probably mostly a result of small number 
statistics; even at high E, there exist regions with no em- 
bedded protostar (downward-pointing triangles in Figure 
[7]). In addition, one can see evidence for correlations be- 
tween N and E for individual regions with reasonable 
stellar statistics (for example, Oph, Cha, and Ser-Aql; 
see Figure[l]for symbol identification). Nevertheless, the 
slope of N/A vs. E is basically determined by the slope 
of the A(E) relation. 

It is clear from the right panel of Figure [7] that a large 
fraction of the scatter in the A(E) relation is due to 
the superposition of different regions (coded by differ- 
ent symbols), each of which having intrinsically similar 
slopes but differing offsets. This leads to the plausible 
conjecture that more massive clouds have larger areas at 
a given surface density, leading to correspondingly more 
stars per E bin. To investigate this possibility, in Figure 
[8] we plot the area divided by a scaled total mass from 
Table 1 of HIO. The result is a greatly reduced scatter 
around the trend line. 

Formally, the least-squares fit using the errors from 
HIO and eliminating the Ser-Aql region for the unsealed 
data in Figure [7] yields 

logA(pc^) — 

(-3.58 ± 0.21) log(E/120Mo pc'^) + 0.56 ± 0.05 . (11) 
while for the scaled data of Figure |8] it is 

logA(pc2)/Mtot(2OOOM0) = 

(-3.24 ± 0.18) log(E/120Mo pc^^) + 0.62 ± 0.05 . (12) 

Here we have chosen to normalize at the approximate 
surface density "threshold" E^ = 12OM0pc"2 of LIO 
and HlOB The power-law behavior ^ E ^ is 
reasonably consistent with that observed in studies of 
the probability de nsity function of column density (e.g., 
iKainulainen et al.l [2 009: Frocbri ch fc RowlesI [2010 ) . al- 
though lFroebrichfc Rowlea (20101) suggest that the slope 
may vary depending upon spatial resolution. 

This analysis of A(E) leads to another important is- 
sue; specifically, the important physical variable is the 
local volume density p, not E. There is no unique way 
to convert a surface density to a volume density distri- 
bution without using density diagnostics. However, it 
is instructive to make the following guess. We suppose 
that the area projected upon the sky A at a given value 
of E is representative of that which would be seen from 
an orthogonal projection, at least on average. In other 
words, we assume that the characteristic average scale 

^ This is not the Larson (1981) scaling relation, p oc 
r~^'^, which would imply S ^ constant (see discussion in 
IBallesteros-Paredes eFaLllSoTl 'l 
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length along the line of sight is h ^ A^/^; thus 

p-S/A^/^ (13) 

Using equation [T^l the value of the critical surface den- 
sity area is A(Ec) ~ 4pc^(Mtot/2OOOM0); this implies a 
median length scale for Mtot ~ 2OOOM0 of /i ~ 2 pc. For 
Ec = I2OM0 pc~^ = 2.5 X 10~^gcm~^, or a column den- 
sity N{H2) — 6.3 X lO^-'^cm"^, the average volume den- 
sity is then n{H2) = lO3cm-3(Mtot/2OOOAf0)-°-5. We 
might go further and use the mass range in Mtot ex- 
plored by HIO to estimate upper and lower ranges. This 
results in estimate densities at Sc of n{H2) ~ 500cm~^ 
and ~ 2.6 x lO'^cm"^ respectively. This is important to 
keep in mind when evaluating models which explain the 
observations of q (§4). 

This analysis of the relation leads to another 

important point; in proceeding to higher values of S, 
HIO is probing smaller volumes. Indeed, applying the fit 
obtained for the median dependence of A{Ti) (equation 
[T2)) . at E ^ IOOOMqpc"^ the typical length scale being 
probed is /i ^ 0.1 pc and thus the implied average vol- 
ume density is n{H2) ~ lO^cm"^, essentially protostellar 
core volumes and densities. 

This interpretation of the linear SFR vs. E relation as 
simply observing gas which is forming stars at maximum 
efficiency is further supported by the fact that the use of 
infrared l uminosity t o det ermine the SFR for the HCN 
clumps bv lWu et al.l (|2010l ) is only possible as reasonably 
massive stars have already formed in these regions. If by 
using HCN a specific density range is being identified, 
and thus a specific free-fall time, a linear relation between 
the star formation rate and gas mass follows. 
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Figure 1. Left: The correlation of the number of protostars per unit area as a function of S, using the data of HIO (their Table 2). 
Differing symbols indicate differing clouds, labeled by the identifiers used by HIO. Downward pointing large open triangles show "upper 
limits", regions where no protostars were found, and so a value of one protostar was assigned for plotting purposes, as in HIO. This is 
essentially the data plotted in Figure 8 of HIO, except that we have not divided by a constant star-forming timescale of 2 Myr, we have 
added the number of Class I and flat-spectrum sources together, and did not include data from the HCN clumps (see text). Right: The 
total protostellar mass (as estimated from the total number of Class I and flat-spectrum sources, multiplied by an average mass of 0.7Mq, 
and divided by the gas mass within the contour, as a function of S. Note that the slope is not significant because the quantities in both 



axes depend upon E (i.e., 
(see text) 



Ma 



S X A(E)). However, it is suggestive that, at high S, the ratio of stellar to gas mass begins to exceed 10% 
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Figure 3. Left: Evolution of mean ring densities for the standard model with Mc = IO'^Mq (solid curve) and simulations with initial 
masses of 1.5 x 10^, 2 x 10"^, 3 x 10"^, 6 x 10"^, 1 x 10'' and 2 x IO'^Mq. The filled points show the time when collapsing fragments appear 
in the ring. Right: evolution of density for the uniform free-falling sphere. 
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Figure 4. The ratio q of gas mass Mg above He to the estimated mass of young stars M» , plotted as a function of the ratio of Mg to the 
total cloud mass, using data from Table 2 of LIO. A mean mass of 0.7Mq per star has been used to convert the number of stars to a total 
stellar mass. Errorbars reflect only counting statistics for the number of stars. There is an indication that q decreases as the fraction of 
dense gas mass increases (see text) 
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Figure 5. Results for the solution of equations [8] and (9] and adopting Mg o = for two values of e. The dotted horizontal line indicates 
the observed upper q, while the two solid horizontal lines indicate the empirical limits for the "evolved clouds" (see text). 




Figure 6. Inferred ratios of stars to dense gas, using the standard model simulation and various assumptions about star formation rates. 
On the left, the evolution of dense gas with a time offset. In the right panel, we show q = Mg/Mt, as a function of time assuming 
Mg = M^e„a(, — Mt , where Mjg„gg is the total gas as in the left-hand panel, and M* the mass in the circular filament above that critical 
density 0.35 Myr earlier, times an efficiency factor e. The dotted curve assumes e = 1, the solid red curve e = 0.3, and the dashed curve 
e = 0.1. The blue solid curve, labeled /, shows the situation assuming that after every 0.35 Myr a new, independent dense gas region forms 
that does exactly the same thing as the e = 0.3 case, just starting later (see text). 
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Figure 7. Left: relation between the number of protostars within a contour of area A as a function of E. The correlation is not strong, 
though individual regions (Ser-Aql, pentagons; Oph, open triangles; Cep, crosses; Per, filled triangles) show a trend of A' increasing with 
S. The weak correlations are probably mostly due to the small number statistics, as evidenced by the contours containing zero protostars 
(downward-pointing open triangles). Right: area of each contour vs. S of that contour. There is a clear correlation; moreover, individual 
regions show similar slopes with vertical displacements. The results indicate that the correlation seen in Figure [T] is driven mainly by the 
A(S) relation. Symbols as in Figures 1 and 8. 
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Figure 8. Areas at differing extinction contour levels, normalized by the scaled total cloud mass above an extinction threshold Ay = 2 
(from Table 1 in HIO) for the clouds in Table 2 of HIO with more than one young star contained within the contour. The dotted line shows 
the least square fit. 



